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Abstract 

co" 

We study a model describing some aspects of the dynamics of biopolymers. The models 
involve either one or two finite chains with a number N of sites that represent the "units" of 
a biophysical system. The mechanical degrees of freedom of these chains are coupled to the 
internal degrees of freedom through position dependent excitation transfer functions. We 
reconsider the case of the one chain model discussed by Mingaleev et al. and present new 
results concerning the soliton sector of this model. We also give new (preliminary) results in 
the two chain model in which case we have introduced an interaction potential inspired by 
C*~) . the Morse potential. 
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1 Introduction 



Recently some work has been done on the study of curvature-induced symmetry breaking 
effects in nonlinear Schrodinger models |1J. The idea here is to study a one-dimensional 
'"O ' discrete Nonlinear Schrodinger equation 

c : 

o . d 

ij^i + Y, J ik^k + x\H 2 *Pi = o (1) 

k 

in which the excitation transfer function (Jik) depends on the positions of the lattice 
points. Often f| one restricts oneself to the nearest neighbour approximation in which 
case Jik = J5i_k,±i but, as pointed out in [I], a more general case involves = J(\fl— fk\) 
where fl describes the spatial position of the i th lattice site. Of course we expect J(\fl— fl\) 
to be a fast decreasing function of its argument. When all the points of the lattice lie 
along a straight line, i.e. fl = a + otib, where a, is a linearly growing function of i, the 
results are not that different from the case of a regular lattice with the nearest neighbour 
approximation. However, as pointed out in pQ , any curvature in r*j can induce extra effects 
which do affect the behaviour of the ip fields. 

In P the authors have studied the case when the lattice points lie on a parabola with 
some curvature. Then they have considered J(a) to be given by J(a) = J exp [—ad) 
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potential whose shape depends on the curvature of the parabola. This has an implication 
on the ground states of the theory. 

This idea has then been applied to the study of the dynamics of biopolymers p j . The 
authors of [Sj have considered the case of an excitation field ip on a one polymer chain. 
The links of the polymers have been allowed to move, and their motion has been controlled 
by their inter-link forces and also by the force due to the excitation field ip. 

The total Hamiltonian describing the system they have considered is given by 

H = T + U + V (2) 
where the kinetic energy T and the inter-particle interaction potential U read: 

T = ¥E(f) 2 and U(f i )=U s + U R + U B . (3) 

i 

The streching energy Us is given by: 

^ = fE(l^-^il- a ) 2 - ( 4 ) 

i 

a is the elastic module of the streching rigidity of the chain and a denotes the equilibrium 
lattice spacing. The repulsive potential 

U R = I E E( rf " 13 " r k \) 2 Q(d - \n -f k \) (5) 

i kj^i 

comes only into play if the distance between two particles becomes smaller than the 
diameter d of the particles themselves. I.e. if the particles start to overlap, they are 
repelled. However, since the strength 5 of the repulsion is chosen finite, particles are still 
allowed to overlap to some extend. The final bit of the inter-particle interaction describes 
the bending energy: 

^ = fE i-( fl X~)' . *i = m- (6) 

i 

k is the elastic module of the bending rigidity while 0{fi) is the angle between the two 
vectors (r*j — r*i_i) and (r i+1 — fty (for a detailed formula see further in this paper or 0). 
6max is the corresponding maximal bending angle. Finally the potential of the complex 
scalar field tp reads: 

v = £{ 2 I^I 2 - £ JftViV* - Ixl^l 4 } (r) 

where x denotes the self-trapping nonlinearity and the exciation-transfer coefficients are 
given by: 

J ik = (e /3 -l)exp(-/5|f i -f fc |). (8) 

Note that the non-linear Schrodinger equation for the if) field and the Newton equations 
for the chain are coupled only by the Jik term. 

In biophysical systems the excitation ipi can be thought of as an amide-I vibration, 
i.e. an excitation in the C=0 bond of the peptide group, which transmits along the 
chain and deforms it. Through this interaction a new localized energy state is created, 
the Davydov soliton 4J, which can transport energy without dispersion. Similarly, the 
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electron-phonon interaction jHE]. 

One end of the polymer has been set to be free; the other one, assumed to be far away, 
has been kept fixed. 

The authors of |3j then have performed many interesting simulations of the equations 
which govern the dynamics of the electron and the lattice. Their most spectacular results 
have involved them showing that a single excitation of the ip field at the free end results, 
after a long period of time, in a gradual folding of the chain. Then they have explained 
their results as coming from the instability generated by the development of the curvature 
of the chain. 

These exciting results have made us think of systems involving more chains with an 
interchain interaction. Proteins often consist of two or more polypeptide chains and take 
on characteristic structures in 3-d space. Examples of such secondary structures are the 
a-helix and the /3-sheet, where in the latter pairs of chains lie side by side being stabilised 
by hydrogen bonds between the carbonyl oxygen atom of one chain and the NH group of 
the other. Furthermore, DNA consists of two polynucleotide chains which are connected 
to each other by hydrogen-bonding and are coiled around each other in a double helix. 

In view of these points, we have decided to extend the investigations of [3] to systems 
of two chains with the simplest interchain interaction as given in, e.g., [0] and inspired 
by the Morse potential [7j. So we have looked at a system of two chains with coupling 
constants given by the appropriate generalisations of those of [Sj. Note that the coupling 
between the ip fields on the two chains allows the field to spread between them altering 
the values of the effective coupling constants on each chain. 

Our model is presented in the next section. 

To test our program, we have, first of all, tried to reproduce the results of |3j. Unfor- 
tunately, ref |3| does not give all the details of their simulations - so in fact, not being sure 
what their initial conditions had been, we have not been able to reproduce their results. 
However, we have found some interesting properties of the system which we discuss in the 
following section. 

Then we have looked at the system involving two chains. Our results are presented in 
Section 4. 

We end the paper with some conclusions and our plans for the further studies. 



2 Our model 

As our model we take a straightforward generalisation of the model of jSj to which we 
have added an interaction between the chains Wi nt . 
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Thus the total Hamiltonian is given by : 

H = T + U + V + W mt (9) 

with the kinetic energy 

^ = EfE(^) 2 > (io) 

3=1 i 

and the inter-particle interaction: 

u(nj) = u s + u B + u R (n) 

where the stretching energy is a sum of the respective stretching energies on chain j = 1 
and chain j = 2: 

Us = jz a ~f Ed*** - ri-i,j\ ~ «i) 2 • (12) 

j=l i 

Similarly, the bending energy reads: 

u B = E f E M JL)- ' % = %(^i) • (1 3 ) 

Finally, the repulsive potential is a sum of the respective repulsive potentials on the two 
chains and a new term, which describes the repulsion if the two chains come closer than 
dj/j to each other: 

2 

U R = Ed E E( rf i - 14' - ^j|) 2 ©(di - ~ r kj \) 

j=l i k^i 

+ Elr 1 Efe - 1^ - v l) 2 )0fe - 1^ - v I) • (1 4 ) 

The energy of the excitation is given by : 

f = EE{2|^f -E4^ - EEWi/ - §x,#d 4 } (is) 

with the energy transfer coefficients: 

4 = ^(e ft - 1) exp(- / S i |f ij - f fc ,|) , (16) 

Klf = (e*r - l)exp(- 7 |7% - . (17) 

The second term on the rhs of (|15J) describes the energy transfer along one chain, while 
the third term corresponds to the energy transfer between the two chains. 

For the interaction potential between the chains we take: 

win* = E E^{ ex p(-«l^- - vl) - 1} 2 - (is) 

i' +3 1 
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arises in a convenient model for the potential in a diatomic molecule. D then is the 
potential energy for the bond formation and a a parameter controlling the width of the 
potential well. Since the minimum of (fT8|) is given for fij = r-' for all i this potential 
leads to an attraction between the chains. 

Note that our terms are simple generalisations of the expressions from [Hj. All fields 
and coupling constants have an extra index j = 1, 2 which tells us which chain they refer 

to. We have also added a term J2i V^Vv coupling tp^ fields on two different 

chains and, as we have mentioned before, Wi nt . In addition, we have multiplied J( k by 
an extra constant A-,- to extend the model to cases where the extension of the interaction 
over the chain and the strength of the interaction can be chosen independently from each 
other. 



2.2 Equations 

It is easy to derive equations which follow from our Hamiltonian. They are given by: 



2.2.1 The Schrodinger equations 



= 2tfj ij - £(e& - 1) exp(-^|r% - f kj \)^ kj 

k+i 

- EE( e7 - 1 ) ex p(-7l^- v I ) %' 

i'+i 1 

- XjlAjl^ij, J = 1,2 (19) 
and, for the chains themselves: 



2.2.2 The Newton equations 



i 

+ EEE^;%-^ > i = w ( 2 o) 

m jVj 1 

where Mj is the mass of the particles on the j-th chain and Uj are the corresponding 
damping parameters. The derivatives of the different potential terms with respect to re- 
read: 
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dUg _ lJ dr tj / 99 x 

drij - K 3 O-(9ij/0 m am) 2 ) 2 \ I 

= -5, J2(dj - Irli - r kj \) ■ ^ - E S fj (d fj - \f, - r^i) ■ ^ (23) 

E E E #C>A-' = 7(1 - O E E «p(-7l^ - v I) • • + ^) 

(25) 

^ = _2Da E{ ex P(-«l^' - Vi) " 1} " S4[ ■ exp(-a|f^ -f^l). (26) 

As in |3j, we have introduced a viscous damping Vj resulting from the assumption of 
having an aqueous environment. This is a physical reason for the introduction of such a 
term. On the other hand we will try to find the lowest energy (stationary) state of our 
system and we will use the absorption to reduce the energy and to 'drive' the system to 
its ground state. Thus, in practice, we will start with a good approximation to the ground 
(lowest energy) state and then evolve it with nonzero absorption - thus letting it 'flow' 
towards this true ground state. 

To proceed further we need to determine the dependence of angles 0^ on fy. Geomet- 
rical reasoning gives: 

9 t3 =n-9 (27) 

where 

9 = arccos( ^^ + | 1 ^^r^;^r 1 - j+T V 1 ; J ^- 1 ' J ). (28) 

Then 

% = + Sp-i + (29) 

The first term on the rhs is given by: 

\d9i£] _ (-1 _ / ^j-^+i,j-^-^-^-i,i+^+i,j-^-i,j \2\-l/2 F^fjj) / qn \ 
Idnjli - K L y [*^-^-wlin+i,i-nj| ' > ' Fa(fa) { } 

F x{ni) = ( 2 nj - r i+1>j - r^^dfij - f^ ltj \ \r i+ljj - fy|) 

■ ( 31 ) 

^2 (ry) = (\rij ~ n-u | |f i+ ij - fijl) 2 . (32) 



with 



and 



The second and the third terms on the rhs of ()29jl are given by similar expressions 
with functions F 3 , F 4 , and F 5 , F 6 , respectively. 
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We have also considered different terms to describe the effects of the curvatures of the 
chains. Thus we have replaced Ub by Uc- 



0c = XXE*S (33) 

j = l i 

where ku is the curvature of the jth chain at site i and so Uc is given by 



2 

" A . X^fl 

\r i j—r i -ij\\r ? i+1 j-r i j\ 



^ = EA J E(l-fctefefif) (34) 



3=1 



where Aj = 2Aj. 



Replacing ^fi- by ^jpf, we get: 

%F- = A,-(M4 + ##4 + #P4) (35) 
with Fi etc. given by the previous equations. 

Now, comparing with the expressions we have had before, we see that 

fcj oc (1 + cosfl) = (1 + cos(tt - %)) = (1 - cos%) (36) 
so that for 9^ = (straight line), fc?- = 0. 

We can also introduce torsion. To do this we note that from Frenet's formulas we have 

r 2 =S-* a =0-« a . (37) 

Thus it would make sense to use an analog of (J35j) and so write 

r 2 oc(l-cos0) (38) 

where <fi is the angle between the two planes we compare. For cf) = it is clear that r 2 = 0. 
In general <fr is given by: 

cos d) = ttfM^ft^jWZi+ld^?^ (39) 

sin0y sme i+1 j\rij-r i - 1 j\\r i+1 j-r i j\ 2 \r i+ 2,j-r i+ ij\ " v / 

This can be rewritten as 



[(nj-n-l,j)-(^+l,j-^j)][(n+l,J-^j)-(n+2,j-n+l,j)] 

sinfly sinfl i+ i j j|fij-Fi_i i j||Fj + i,j-f ? ij| 2 |ri + 2,j-ri + ij| 

[(rij-ri_i,j)-(r i+ 2,j-ri + i,j)] 
sinfly sin0i + ij|r* i:) -ri_i J Hri + 2,j-ri + i ij | " 



(40) 



We know that sin 6^- = sin 6^- and that sin(arccos(a;)) = yl — x 2 . Thus 



= v 1 - ( v ^-g^;i j,r - ij ) 2 ( 4i ) 

For a straight line, we have to be careful since in Eq. (j39|) sin^- = 0. 
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and of the excitation field. The preliminary results indicate that the dependence on the 
detailed form of these terms is not very strong. So in this paper we report the results for 
the case of Ub leaving the detailed study of the dependence on the form of Uq to a future 
publication. The same applies to the effects associated with non-vanishing torsion. At 
the same time our observation of the weak dependence on the details of Uc shows that 
our results are quite generic in their nature with most observed effects determined by the 
other terms in our Hamiltonian (2). 



2.3 Initial and boundary conditions 



Next we have to decide on the initial and boundary conditions. 

First of all we have chosen the chains, initially, to lie parallel to each other (and we 
have placed them along the x axis in the xy plane, starting at x = and running to x = N, 
where N is the total number of links of each chain). So we have put r iX \ t=0 = (xn, 0, 0), 
and r i2 \t=o = (2^2 , 1,0). It can then be shown that all angles vanish identically: 

9\ t =o = arccos(-l) = vr -> %| t=0 = 0. (42) 



It is possible to fix the ends of the chains (x = N) by requiring that 

dt u ' dt u ' dt u - y^°) 

However, in most of our simulations, the end has been allowed to move freely. 

We cannot define 9y at x — 0. Thus, at this point, we have taken the bending angle 
as the angle between the x-direction and the vector r 2 j — fy : 



y = arcsin( v / (y 2j - y^) 2 + (z 2j - z^f/\r 2j - r y |). (44) 



3 One chain 



In this case we have put cr 2 = k 2 = 82 = ^12 = D = j3 2 = 7 = X2 = ^2 = 0. Thus the 
second chain is completely decoupled and our problem is reduced to that of ref [3]. 

We have redone some of the work reported in |3j . We have started the simulations by 
putting, initially, ifj^i = 1 and ipi,i — 0, when i = 2,..N. Using the parameters of 
ie a x = 1000, k x = 0.06, 8 X = 100, (3 X = 2, X i = 3.2, v x = 0.3, 6 max = vr/3, M x = 0.5, 
d x = 0.6 and in addition Ai = 1, we have found that the links of the chain have moved 
very little, and that the motion was only in the x direction. This is easy to understand 
as nothing breaks the symmetry keeping the motion restricted to the y = z = line. Of 
course, this symmetry can be broken dynamically through the numerical inaccuracies but 
this process is extremely slow. 
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merit of the chain in the y as well as in the z direction. Thus we put initially y^i = 0.1, 
zi t i = 0.01 and set all other y^i, z it i equal to zero. In all our simulations, the energy has 
decreased due to the absorption and has settled quite quickly to its final value. 

We have first studied the existence of a soliton-like structure as a function of the 
parameters (3\ = (3 and Ai. We have found that the range of Ai for which the soliton 
doesn't get destroyed is very limited for all (3\. For Ai > 5, ijj^i is completely spread over 
the chain after some finite time, while for Ai = 2,3, a soliton-like structure is still seen to 
move up and down the chain, however, a lot of "background" noise is present. A sharp 
soliton seems to be present only for Ai ~ 1. We have thus fixed A x = 1 and studied the 
dependence of the height, width and velocity of the soliton on the parameter (3\ = (3. For 
the position i max and the height (ipi,iip*i) m ax = l^ilmax °f the soliton's maximum, we 
have used a quadratic approximation. For the width A of the soliton in the x-direction 
we can use two different definitions: 

• a "quantum-mechanical" definition: 

A qm = V< x 2 > - < x > 2 (45) 

with < x 2 >= ^Ei^,il^,i| 2 and < x > 2 = £j x;,i|V>i,i| 2 ) 2 > or 

• a definition using the half- maximum of ijjijj*: 

A V2 = x + 2 _ x - 2 (46) 

where x\, % > x max and Xy 2 < x max denote the two x-values for which \i>i,i\ 2 (xi/ 2 ) — 

IV^I | 2 (zi/2) = (\^i,l\ 2 )max/^- 

Our results for the height, position and width of the soliton as a function of /3x = after 
t = 26.4 sec are shown in Fig.s 1(a), (b), (c), (d) for three different values of the coupling 
Xi = X- From these figures we see that the soliton moves quicker with the decreasing 
of (3. Moreover, the height of the soliton's maximum decreases with the decreasing of (3, 
while at the same time the soliton gets broader as can be seen from Fig.s 1(c) and 1(d). 
Comparing Fig.s 1(c) and 1(d) it is also clear that the qualitative results from the two 
different expressions for the width are in good agreement with each other, even though the 
"quantum mechanical" expression takes into account all tpn for i < 40 along the chain. 

For (3 > 10, the excitation is completely trapped at its initial position = 1. This 
can be explained by noting that for large (3, the ip field equation essentially decouples from 
the chain positions. For < 1 it spreads over the chain with an oscillating behaviour of 
l^il ■ This is already noticeable for (3s slightly larger than this value, where the soliton 
maximum decreases as it moves along the chain. For the reasons above, we can only 
calculate the maxima of the soliton etc. for the interval (3 e [1 : 10]. 

It had already been noticed in |H] that the soliton exists only for 2.5 < x < 3.5. We 
find the same in our simulations with the initial excitation being trapped at its initial 
position for x > 3.5 and being completely destroyed for x < 2.5. Thus we have chosen 
X = 2.8, = 3.0 and = 3.2 to investigate the dependence on this parameter. Clearly, as seen 
from Fig. 1(a), the soliton moves faster as x decreases. Moreover, its height decreases 
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for x = 3.0 is not very conclusive; however, we believe that this is due to the quadratic 
approximation we have used. Both methods of computing the width of the soliton again 
give the same qualitative result that the width decreases with increasing x- Thus the 
soliton is "sharper" for larger values of x as could have been expected. 

Following [B], we have further studied the correlation between the position of the 
soliton and the average displacement of the site i given by : 

S i+ i - Si-x = \jdx 2 i+1)1 + dyf +l)1 + dz? +1>1 - yjdx?_ ltl + dy}_ xl + dzf_ 1A (47) 

where dx^i, dy^ and dz^i are the displacements from the initial position of the site i in 
the x, y and z directions. For two typical examples with x — 3.2, /3 = 5 and (3 = 10, 
respectively, with all other values given as before, we show Si as well as 1^,1 1 2 i n Fig.s 
2 (a) and 2 (b). Clearly there is a correlation between these two quantities in the sense 
that the maximum of | 2 is located in the region where Si has a local minimum. Of 
course, the soliton movement generates also some "background" noise which influences 
the displacement of the sites from their initial position. The minimum of the displacement 
at the left end of the chain can be attributed to the fact that we have introduced an initial 
displacement of the chain in both the y and z directions. 



4 Two chains 

We have performed several studies of a two chain system. In each case we have originally 
placed two chains along the x > axis, one of them at y = 0, the other at y = 1 (with 
z = 0). Such a configuration does not depend on z and this symmetry is preserved by 
the equations of motion. Hence, to see any effects like folding, we have chosen to break 
the symmetry explicitly by displacing, initially, the first two links of the first chain in 
the z direction by: z\,\ = 0.2 and z 2 ^ = 0.1. We have also used other values for this 
displacement but we have found our results not to be too sensitive to the initial values of 
this displacement (as long as it was nonzero). 

Then we have considered several initial conditions for the ipij fields. In particular 
we have performed many simulations with ■012,1 — 1 (and all others are zero), or with 
■012,1 = 0.95, ip!3 t ! = i0.3 and all other vanishing. 

We have performed these simulations for many values of f3j and \y F° r the other 
parameters we have used ox = a 2 = 3000, a x = a 2 = «i = k 2 = 1, S\ = S 2 = S\ 2 = 100, 
d\ = d 2 = d 12 = 0.6 and 7 = 2, d — 0.5, a = 1.8. 

At first this evolution is very fast, then it slows down and after a while it approaches 
the final configuration exponentially slowly. During the evolution at first the changes of 
the positions of the chains are clearly visible, then become hardly noticeable and finally 
the chains look as if they were essentially static. The ip fields still evolve but even their 
evolution is not very dramatic. Having reached this stage we are reasonably confident 
that our configurations are not very different from the final asymptotic lowest energy 
configuration. In fact, in each case, to be absolutely certain that we are well "beyond" 
any transient effects, we have run each of our simulations over periods of several months 
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results. They are somewhat qualitative. 

We have preformed several simulations (for a range of values of Xj (from 2 to 9) and 
for various values of f3j (from 1.5 to 2.5). 

All the simulations have shown a localised (soliton-like) behaviour of the ipij fields, i.e. 
the effects on ipij fields were reasonably well localised. This localisation concerns both 
chains and it depends crucially on the values of Xj- F° r larger values (say Xj — 6.4) we 
end up with a well defined soliton-like peak (clearly visible when you add the two chains) 
oscillating (relatively little) in size. Like in the case of a single chain the soliton moves 
and its speed depends on f3j. For larger values of Xj the soliton is narrower and so more 
sharply defined while for smaller values it is more spread out. For Xj — 3.2 one can still 
identify a clearly localised structure but whether this structure should be called a soliton 
remains to be decided. 

We have found that the dynamics of the chains themselves depends strongly on the 
soliton structure. For the case of small Xj an d thus spread out \ipij\ 2 , the chains get very 
much deformed (at least over short periods of time) in the sense that they have large local 
curvature. For sharp solitons however, i.e. large Xji the behaviour of the chains is very 
similar to that observed in the one chain model. Namely, the chains change their shape 
very little and remain more or less straight. 

We show a typical example of the minimal energy configuration in Fig. 3. This is for 
fli = @2 — 2.0 and xi — X2 = 4.0. Note that this example corresponds to a spread out 
l^ijj 2 . Clearly, the chains have large local curvature for < x < 60. Looking at one 
chain by itself, it seems that loops have formed. Especially the red-coloured chain has 
formed a large one. Moreover, at small x, the two chains have overlapped each other in 
the x-y-p\ane and it seems like they are winding around each other. This winding also 
happens in the a-helix, the feature that appears e.g. in DNA. 

Thus our results suggest that solitons exist when the chains are little deformed; while 
the large deformations of the chains correspond to more spread out \ipij\ 2 fields. Moreover, 
we have not observed any behaviour which suggests any "folding up of the chains" . At this 
stage our results are a little qualitative. We hope to present a more quantative discussion 
of our results and an analysis of the dependence on f3 and k in the next paper (sometime 
in the future) [Oj. 



5 Conclusions and Outlook 

In this paper we report our first results on the dynamics of biopolymer chains. While the 
authors of [3] put emphasis on the dynamics of a chain itself (they studied one chain) and 
found that for suitable choices of the coupling constants the chain folds up, we have been 
mainly intersted in the dynamics of the soliton moving along the chain. We have found 
that the soliton exists only for very specific choices of the coupling constants. We have 
studied numerically how the speed, height and width of the soliton depend on the coupling 
constants. Morever, we have confirmed the existence of a direct correlation between the 
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For the case of the two chain model, we have introduced suitable generalisations of 
the energy transfer coefficients and the rigidity potential. Moreover, we have added an 
attractive interaction potential between the chains inspired by the Morse potential. We 
have found again that the existence of the soliton depends crucially on the non-linearity 
parameter x- However, while for one chain the soliton is either completely spread over 
the chain or sharpely localised, it seems that in the case of two chains there exist a 
sort of intermediate situation. In this, a number of small peaks in \%p\ 2 exist that travel 
up and down the chain. Remarkably, the dynamics of the chains themselves becomes 
only interesting in this latter case. We find features like the formation of loops and the 
winding of the two chains around each other which also happens in biological systems like 
e.g. DNA. We plan to report a more detailed quantative analysis of the two-chain model 
in a future publication [Hj. 
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Position of soliton after t=26.4 sec 
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Figure la: The dependence of the position of the soliton's maximum at site i is shown as 
function of (3± = f3 after t = 26.4 sec for three different values of %• 
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Height of soliton after t=26.4 sec 




Figure lb: The height of the soliton's maximum (ipil>*)max after t = 26.4 sec is shown as function 
of Pi = (3 for three different values of %. 
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Width A 1/2 of soliton after t=26.4 sec 




Figure lc: The "half-maximum" width A 1 / 2 of the soliton after t = 26.4 sec is given as function 
of Pi = (3 for three different values of %• 



15 



Width A qm of soliton after t=26.4 sec 
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Figure Id: The "quantum-mechanical" width A qm of the soliton after t = 26.4 sec is given as 
function of f}\ = f3 for three different values of %• 
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Figure 3: The minimal energy configuration of the two chains for (3i = (3% = 2.0 and Xi = X2 = 
4.0 is shown in x-y-z-space. The two dashed lines at z = —0.5 represent the x-y-plots [z = 0) 
of the chains. 
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